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Abstract: Cosmological magnetic fields are being observed with ever increasing corre- 
lation lengths, possibly reaching the size of superclusters, therefore disfavouring the con- 
ventional picture of generation through primordial seeds later amplified by galaxy-bound 
dynamo mechanisms. In this paper we put forward a fundamentally different approach that 
links such large-scale magnetic fields to the cosmological vacuum energy. In our scenario 
the dark energy is due to the Veneziano ghost (which solves the U(1)a problem in QCD). 
The Veneziano ghost couples through the triangle anomaly to the electromagnetic field 
with a constant which is unambiguously fixed in the standard model. While this interac- 
tion does not produce any physical effects in Minkowski space, it triggers the generation 
of a magnetic field in an expanding universe at every epoch. The induced energy of the 
magnetic field is thus proportional to cosmological vacuum energy: pem — B 2 ~ (£^) 2 Pde, 
Pde hence acting as a source for the magnetic energy pem- The corresponding numerical 
estimate leads to a magnitude in the nG range. There are two unique and distinctive 
predictions of our proposal: an uninterrupted active generation of Hubble size correlated 
magnetic fields throughout the evolution of the universe; the presence of parity violation 
on the enormous scales 1/H, which apparently has been already observed in CMB. These 
predictions are entirely rooted into the standard model of particle physics. 
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1. Introduction 

The origin of cosmological magnetic fields, that is, magnetic fields which permeate the 
largest structures found in the universe such as galaxies, clusters, and so on, still enjoys 
the status of "yet to be solved" mistery. Such cosmologically correlated magnetic fields have 
historically been discovered at ever increasing distances (see the comprehensive reviews [1, 

2, 3, 4, 5, 6], mostly due to the fact that experimentalists have come up with new methods 
and new devices capable of observing their effects in the ever farer and bigger structures 
in universe. 

The models and mechanisms that aim at explaining the origin of the large-scale mag- 
netic fields can be roughly filed under two wide categories, astrophysical and cosmologi- 
cal ones. The astrophysical mechanisms mostly rely on what is typically termed battery, 
see [2, 6] for a review. Such mechanisms generally provide protogalactic magnetic fields 
but they are not likely to be correlated much beyond galactic sizes. 

On the other hand, cosmological models for the origin of cosmological magnetic fields 
typically involve some sort of (possibly first order) phase transition, such as the electroweak 
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transition or the QCD confinement one, in order to generate Hubble scale seed fields which 
then are amplified through the mechanism known as dynamo [5, 6]. The typical correlation 
length of such seeds is typically many orders of magnitude smaller than needed today, unless 
a process called inverse cascade takes place (see for instance [7, 8, 9]). An inverse cascade, 
in brief, spreads highly energetic short wavelength modes to larger correlation lengths at 
smaller amplitudes. This effect is expected to take place whenever a helical field evolves in a 
turbulent plasma, as the early universe mostly is, see again [2, 5, 6] for details. Depending 
on the architectures used to model such phenomena, one expects to obtain correlation 
lengths that vary between a few pc to the very optimistic 100 kpc, but a very weak field 
intensity. 

In the case of an inflationary universe there is also the opportunity of a super-horizon 
magnetogenesis due to the peculiar expansion character of an inflating universe [10], which 
however demands a beyond the standard model (SM) coupling for electromagnetism for 
the latter is Weyl invariant and is not amplified by gravitational interactions (unless one 
resorts to the conformal anomaly, see [ll]-but that again is inefficient as long as we stick 
to SM physics). 

In all cases, the dynamo is a necessary step in order to revitalise the weakened magnetic 
fields. The dynamo must be very efficient for most of the models pushed forward, and, more 
importantly, should be operating at very large-scales (beyond galaxies). All these scenarios 
however incur in a series of difficulties when faced with observations; these problems become 
profoundly more difficult as observations reach further in intensity and distance, as we argue 
below. 

The structure of our presentation is as follows. In the following section 2 we discuss 
how observations compare with most theoretical paradigms on the generation of large scale 
magnetic fields. We elaborate on a number of problems which appear to be beyond the 
abilities of such models. Moving on to section 3 we offer an alternative approach to resolve 
the mystery of large scale magnetic field based on the anomalous coupling between dark 
energy and magnetism. In subsection 3.1 we review the dynamics of the Veneziano ghost 
(which solves a fundamental problem in QCD and plays the role of source for the DE 
in our framework) in an expanding universe, and sketch how the observed vacuum energy 
arises, and how the auxiliary conditions on the physical Hilbert space (similar to the Gupta- 
Bleuler condition in QED) keep the theory unitary. In subsection 3.2 we derive the coupling 
between the Veneziano ghost and electromagnetism via the conventional triangle anomaly. 
We shall argue that this coupling, despite being of the same order of magnitude of <7 7r o 77 
(which describes the tt° — > 77 decay) , still does not lead to any physical effects in Minkowski 
space. By contrast, in an expanding universe when the Veneziano sources the appearance 
of the cosmological vacuum energy, it does generate large scale magnetic fields which we 
could be observing now. In section 4 we outline the arguments leading to our quantitative 
estimates on such fields, with particular emphasis on the radical differences between the 
approach pursued here and the far more common generation through primordial seeds. 
Lastly, the final section 5 wraps up the main ideas of the paper with a concise summary 
with conclusions. 
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2. Cosmological Magnetic Fields: Observations vs Theoretical Models 

2.1 Observations: ever increasing correlation lengths 

Large-scale magnetic fields have been first discovered in our Milky Way, and subsequently 
in a number of other galaxies of different sizes and shapes, with characteristic intensity of 
around a few ^G. The correlation lengths of such fields range between a kpc and 30 kpc. 
However, this is not the end of the story, as magnetic fields of very similar strengths have 
been observed in clusters of galaxies, where they appear to be correlated over larger dis- 
tances reaching the Mpc region. It is important to notice that such fields are not associated 
with individual galaxies, as they are observed in the intergalactic medium as well [12, 13] 
(see also the recent papers [14, 15]). This poses a serious challenge to astrophysically based 
mechanisms, for then for some sort of conspiracy a series of unrelated galaxy-bound fields 
would all choose to align in one specific direction. 

As for cosmological magnetogenesis, already at this point one may notice that even the 
most optimistic and efficient inverse cascades are not going to be able to stretch primordial 
seeds to such wide lengths, as far as the seeds come from the electroweak phase transition, 
or even in the case of the QCD one. In the latter case the best result one can hope 
to obtain is around 30 kpc (see for instance [16, 17]) using Son's estimates for inverse 
cascade parameters [9], which are actually already known to be over-optimistic, see [18]. 
Moreover, recent analyses of the development of magnetic fields in cosmological turbulent 
plasma show that more realistic parameters for the inverse cascade lead to even smaller 
correlations [19, 20]. 

The most recent hints towards a possible magnetisation of gigantic supercluster struc- 
tures [21, 22], whose size can easily be two orders of magnitude beyond the clusters within 
it, although not yet fully conclusive, would represent a further theoretical challenge in 
modelling magnetogenesis, pushing the correlation lengths further away up to fractions of 
Gpc. It is clear that in this case all of the mechanisms mentioned above have very little 
chance of being successful: such correlation lengths are simply beyond their capabilities. 

2.2 Observations: high redshifts 

In conjunction with the observations of supercluster magnetic fields there are new pieces of 
evidence that show how fj,G magnetic fields were present at much earlier epochs than pre- 
viously thought [23, 24]. Given this empirical picture, one immediately sees the problem in 
all mechanism which require a very efficient dynamo amplification be brought to successful 
completion, for if strong magnetic fields are observed at redshifts as high as z ~ 5 there is 
no time for the galaxies to perform enough rotations and pump up the field as hoped. In 
addition to that, it looks at least very unlikely that all galaxies amplify initial, commonly 
intense seeds just in the same way to return fiG fields everywhere, for that implies roughly 
the same effective number of turns for all of them. 

2.3 Theoretical models: cosmological seeds? 

Cosmological seeds suffer from another problem besides those just mentioned, as it has 
recently been emphasised in [19]. This can be schematised as follows. Primordial magnetic 
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fields trigger the appearance of gravitational waves which would be around during the 
formation of the light elements (Big Bang Nucleosynthesis - BBN). Such component is 
tightly constrained by the beautiful concordance between theory and observations (a review 
can be found in [25]), and puts problematic limits on the maximal intensity of primordial 
magnetic fields. In particular, for non-helical seeds it has been show that neither inflation- 
born, nor electroweak-born, and not even QCD-born ones would be able to satisfy such 
limits and provide enough power for a dynamo [26, 27]. In the helical case these limits 
are mitigated, but still on the verge of being excluded, for they demand the most efficient 
amplification [19]. 

One further issue plaguing inflation-dawned magnetic gemmae is the fact that the 
back-reaction of the field on the inflaton itself is, in many practical realisations, of non- 
negligible importance [28]. Indeed, such back-reaction would make inflation unable to 
reproduce the observed spectrum of temperature fluctuations as observed in the cosmic 
microwave background (CMB), unless the seeds are smaller that what is required to feed 
galactic dynamos [28]. 

Finally, if the dynamo amplification of primordial seeds serves as the main field gener- 
ation mechanism, it is hard to imagine that low mass dwarf galaxies could produce again 
just the same /iG fields, because of their relative lower rotation rates. 

2.4 Few more hints on existence of large scale fiG fields 

Briefly, let us list a few scattered ideas supporting an ever existing and all-pervading fairly 
strong magnetic field. 

Coles [29] and Kim et al. [30] have argued that sufficiently strong magnetic fields should 
be around during most of the structure formation epoch of the universe, for they would help 
solving some of the discrepancies between ACDM simulations and the observed structures. 

The very simple observation that magnetic fields of intensity spanning at most one 
order of magnitude in the microgauss are found in completely different environments, from 
galaxies to intracluster medium, seems to be telling us that such fields are allegedly inde- 
pendent on the matter density they are found to be immersed in. Indeed, Kronberg [1], 
corroborating this suggestion with a number of other pieces of evidence, had proposed that 
galaxies and clusters have formed in a (iG environmental field, rather than the other way 
around. 

Let us close this overview by mentioning the few models which account for magnetic 
fields in the formation of very peculiar structures such as filaments, which work better than 
their magnetic-less counterparts [31, 32, 33], and the correlation between star-burst rates 
and the intensity of a pervading strong magnetic field within the galaxy, which holds to a 
good accuracy and may signal how stars are given birth in strongly magnetised surround- 
ings [34]. 

3. Coupling Dark Energy with Magnetism 

On a completely different side, naively unrelated to large scale magnetic fields, the last 
decade or so has seen the scientific community steadily realise and acknowledge the exis- 



- 4 - 



tence of the so-called dark energy component which appears to be accounting for over 70% 
of the total energy density of the universe [35, 36, 37]. In this unclustered energy density 
the entire cosmos is immersed, and its origin is very often thought of as the most intricate 
problem of modern cosmology and particle physics. 

In this work, we propose that the two problems are intertwined at their hearts, and 
suggest that the same mechanism which explains the fundamental origin of the cosmologi- 
cal vacuum energy should also be able to encompass the generation of cosmologically sized 
magnetic fields with intensity in the /iG range. More specifically, it has been recently been 
shown [38, 39] how dark energy can be explained entirely within the SM, more precisely, 
within QCD, without the need for any new field or symmetry. In the proposal [38, 39], the 
information about the vacuum dark energy is carried by the so-called Veneziano ghost [40], 
whose properties in the expanding universe will be reviewed shortly. The basic idea is that 
the Veneziano ghost, which is a known, unphysical, degree of freedom in QCD, gives rise 
non a non-zero vacuum energy ("ghost condensation" similar to "gluon condensation" in 
QCD) if the universe is expanding. This effect of the Veneziano ghost is in many respects 
akin to the known phenomenon of particle emission in a time-dependent gravitational back- 
ground, with the important difference that there is no emission of "real", i.e., asymptotic 
states here 1 . Rather, the effect should be interpreted as the injection of extra energy (in 
comparison with Minkowski space) into ghost waves when the universe slowly expands. 
The average momentum resulting from this pumping is obviously zero as momentum is 
still good quantum number in the expanding universe. Overall, this effect is clearly very 
tiny as it is proportional to H. What is also important is that the extra energy is stored 
in the form of ghost waves with momenta ook — k ~ H, as only this much energy can be 
lent by the expanding background, higher frequencies being exponentially suppressed [38]. 
The arguments presented above imply that the typical wavelengths associated with this 
energy density are of the order of the Hubble parameter, ~ \ jk ~ 1/H ~ 10 Gyr, and 
the corresponding excitations do not clump on scales smaller than this, in contrast with 
all other types of matter. 

The very same ghost field couples via the triangle anomaly to electromagnetism with 
a constant which is unambiguously fixed in the SM. While the interaction of the ghost 
with electromagnetic field is sufficiently strong, it does not produce any physical effects in 
Minkowski space as a result of the auxiliary conditions on the physical Hilbert space (similar 
to the Gupta-Bleuler condition in QED) that are necessary to keep the theory unitary. 
However, the induced extra time-dependent energy due to the Veneziano ghost in the 
expanding universe automatically leads to the generation of the physical electromagnetic 
field. What is important is that the typical momentum k EM of the generated EM field will 
be of the same order of magnitude of a typical momentum of the the ghost k. Consequently, 
typical frequencies of the generated EM field will also be the same order, oj em — k EM ~ H. 

One should remark here that our preference in using the approach of Veneziano in describing the QCD- 
related vacuum energy is a matter of convenience. In principle, the same physics is also (hidden) in Witten's 
approach to the resolution of the (7(1) problem [41]. However, the corresponding technique is much more 
involved when the system is promoted to time dependent backgrounds, see subsection 3.1 below, and [38] 
for a detailed discussion. 
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3.1 Dark energy and the Veneziano ghost 

It has been suggested recently that the solution of cosmological vacuum energy puzzle may 
not require any new field beyond the SM [38, 39]. The idea is based on the philosophy that 
gravitation can not be a truly fundamental interaction, but rather it must be considered 
as a low-energy effective quantum field theory (QFT) [42]. The first application of this 
proposal was the computation of the cosmological constant in a spacetime with non-trivial 
topological structure when a Casimir type effect emerges. It was shown that the cosmo- 
logical constant does not vanish if our universe is enclosed in a large but finite manifold 
with typical size L ~ 1/H, where H is the Hubble parameter. The cosmological vacuum 
energy density p\ in this framework is expressed in terms of QCD parameters for Nf = 2 
light flavours as follows [38, 39]: 

2^1^)1 ^ l 4 

This estimate should be compared with the observational value pa ~ (2-3 • 10 -3 eV) 4 [37]. 
The deviation of the cosmological constant from zero is entirely due to the large but finite 
size L of the manifold, and, as we have anticipated, should be understood as a Casimir 
effect in QCD . This proposal has a very simple and analytically trackable analogue in the 
2d Schwinger model [43], and could be tested in upcoming CMB experiments [44]. 

The result (3.1) is a direct consequence of the existence of a very special degree of 
freedom in QCD, that is, the Veneziano ghost [40]. This field is unphysical in Minkowski 
space (it belongs to the unphysical projection of the Hilbert space) in a way akin to the two 
unphysical polarisations of the electromagnetic four-potential in conventional QED. The 
effective Lagrangian for this field was known already in the '80 [45, 46], but it has been 
recently reworked in a very convenient form for studying its curved spacetime properties 
as [38] 



+N f m q \(qq)\ cos 



+ <t>2 



fv' 



(3.2) 



where the covariant derivative is defined as = d^ + T^ so that, for instance D^V" = 
d^V u + r^ A l/ A . The fields appearing in this Lagrangian are 

4> = physical r\ , 4>\ = ghost , 02 = its partner . (3-3) 

It is important to realise that the ghost field 4>\ is always paired up with 02 in each and 
every gauge invariant matrix element, as explained in [38]. The condition that enforces 
this statement is the Gupta-Bleuler-like condition on the physical Hilbert space 7i p hys for 
confined QCD, and reads like 

(02 - 0i) (+) |H phys > = 0, (3.4) 

where the (+) stands for the positive frequency Fourier components of the quantised fields. 
In Minkowski space one can ignore the unphysical ghost field 0i and its partner 02 in 
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computing all S matrix elements precisely in the same way as one always ignores the two 
unphysical photon polarisations when the Gupta-Bleuler condition in QED are imposed. 

However, the requirement (3.4) could not be globally satisfied in a general background 
as explained in details in [38]. This is due to the fact that the Poincare group is no longer 
a symmetry of a general curved spacetime (including the FLRW universe) and, therefore, 
it would be not possible to separate positive frequency modes from negative frequency 
ones in the entire spacetime, in contrast with what happens in Minkowski space where the 
vector d/dt is a constant a Killing vector, orthogonal to the t = const hypersurface, and 
the corresponding eigenmodes are eigenfunctions of this Killing vector. The Minkowski 
separation is maintained throughout the whole space as a consequence of Poincare invari- 
ance. Therefore, all physical effects related to the ghost dynamics are proportional to the 
deviation from Minkowski spacetime geometry, i.e., to the rate of expansion H. This is 
the very reason for the Veneziano ghost to exhibit non-trivial dynamics in an expanding 
universe; gravitational interaction however intervene and change this picture, allowing for 
the appearance of a non-trivial energy density in the time-dependent background. We refer 
the interested reader to the original paper [38] for the details. 

One more comment concerns the appearance of the scale L ~ 1/H in the energy 
density (3.1). As it has been extensively explained in the previous letter [39], and also in 
the longer paper [38], the potential felt by the ghost and its partner is a Casimir-like energy 
which is a result of a subtraction procedure that compares the values of the vacuum energy 
in Minkowski space with that in a general compact manifold (such as a torus of size L) . This 
prescription aims at extracting the physical and measurable portion of the vacuum potential 
energy of the ghost field, by taking such difference between the vacuum energy in compact 
curved space and that in infinite Minkowski space. Essentially, this is our definition of the 
vacuum energy when the "renormalised energy density" is proportional to the departure 
from Minkowski spacetime geometry and remains finite. The basic motivation for this 
definition is the observation (3.4) that in Minkowski infinite space-time the energy due to 
the ghost identically vanishes. Technically, it implies that the Lagrangian itself (3.2) does 
not have any small parameters such as 1/L. However, the vacuum energy thereby defined 
exhibits a Casimir-like effect. Notice that the correction (which was computed exploiting 
the topological susceptibility of QCD when the ghost is present) is linear in the inverse 
size of the manifold, not exponentially suppressed, as one would normally expect in the 
confined phase of QCD where all physical degrees of freedom are massive. 

The ghost we are and will be working with here, and whose effects are central for our 
discussion, was postulated by Veneziano in the context of the U(l) problem. However, 
the same problem had been tackled from a different perspective (although in the same 
large- N c context) by Witten in [41]. In his approach the ghost field does not ever enter 
the system, and makes us wonder whether its consequences are physical or just artificious. 
Without going into the details (see [38] for a technical explanation) , it will be enough here 
to mention that the curved space effects we have directly computed with the help of, and 
attributable to, the ghost, are not going to disappear if we follow Witten's approach. The 
relevant physics will be hidden in the contact term which will depend in a highly non-trivial 
way on the properties of the spacetime (such as curvature) once the apt renormalisation 
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procedure in the expanding background is performed 2 . 
3.2 The ghost-photon anomalous coupling 

Our next step is to include the EM field into the low energy effective Lagrangian (3.2). 
First, in order to do so, we need to know the interaction of the rj 1 field with photons. After 
that, we can recover the interaction of the ghost field <f>i and its companion <p2 with the 
electromagnetic fields because and 4>2 always accompany the physical 7/ in a unique way 
similar to the interaction term (3.2). 

The interaction of the rj' field with photons is a textbook example crafted on the almost 
identical well known anomalous term describing the 7r° — > 77 decay, see, e.g., [47]. The 
only difference is that 7/ is an isotopical singlet state while 7r° is an isotopical triplet. The 
interaction term is 



where a is the fine-structure constant, f^i is the decay constant for the 7/, N c is the 
number of colours, and the QiS are the light quarks electric charges. Finally, F^ v is the 
usual electromagnetic field strength (in curved space), and = e^poF^ /2 its dual. 
We choose e^ vpcT = e^,/ P<J /V — 9 w hh the Minkowski antisymmetric tensor following from 
e M 23 = an d 9 = det g^y the determinant of the metric tensor. 

One should remark here that such kind of anomalous interaction has been studied 
in great details in particle physics as well as in cosmology. In particular, the axion (or 
any other pseudoscalar particle) has exactly the same structure and has been thoroughly 
analysed in cosmological contexts, see e.g., [48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58]. 

We are not really interested in 7/ physics as the heavy 7/ meson of course is not 
excited in our universe. We are interested in the interaction of the ghost field (f)\ and 
its companion <p2 with the electromagnetic fields. The corresponding terms have never 
been discussed in the literature because they do not appear in any gauge invariant matrix 
element in Minkowski space as a consequence of the auxiliary condition (3.4). In fact, 
these unphysical degrees of freedom can be completely integrated out in that case, such 
that they even disappear at the Lagrangian level (this is exactly the procedure that was 
adopted in the original paper [40], see also [38] in the context of the present work). In 
curved space these fields can not be swept under the carpet as they carry relevant physical 
consequences: we must explicitly deal with them. 

In order to derive the interaction term between the fields <pi, 4>2 and one should 
repeat the steps described in [38]. We can explicitly check that the physical rj' always 
enters the Lagrangian in the combination (?/ + 4>2 — <^i)//r?'- Hence, the interaction term 
we are interested in has the structure 



2 The required procedure, for the non-abelian and strongly interacting theory under consideration is not 
known yet. 
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(3.6) 
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For our future discussions we safely neglect the massive physical r( field, and keep only the 
ghost field <f>\ and its companion fa along with the EM field, 

C = -\f^F» u + ^D^D^fc - ^D^D^x (3.7) 

(^)^ + ^*>i«(^). 

where the electric and magnetic fields are the usual Minkowski ones (not rescaled by the 
scale factor of the universe a(t)). We claim that the expression (3.7) is the exact low 
energy Lagrangian describing the interaction of the ghost field <f>\ and its companion 02 
with electromagnetism in the gravitational background defined by the Tu- In Minkowski 
space the expectation value 

<^ph y s|(02-0i)|H phys > = 0, (3.8) 

vanishes, implying that (pi and <f>2 are decoupled from QED, as they should in order to 
preserve the unitarity of the system. 

However, as we have pointed out before [38], the constraint (3.4) can not be globally 
maintained in the entire space in a general curved background. Thus, the ghost field and 
its partner do couple to electromagnetism, and consequently we do expect some physical 
effects to occur as a result of this interaction. Notice that we are not claiming that the ghost 
field becomes a propagating degree of freedom, or becomes an asymptotic state. Rather, 
we propose the description in terms of the ghost as it is a convenient way to account for 
the physics hidden in the non-trivial boundary conditions, see [38] once again. In the 
context of this discussion, a very illuminating example is that of 2d Rindler spacetime 
and the associated Unruh effect, see [59], where all different approaches are workable and 
comparable to give the same result. One should emphasise that the Veneziano ghost is 
very different from all other types of ghosts, including the conventional Fadeev-Popov 
fields. The peculiar feature of the Veneziano ghost resides in its close connection with 
the topological properties of the theory, and expresses the necessity to sum over different 
topological sectors in QCD [59]. This uniqueness manifests itself, in particular, in the 
spectrum of the Veneziano ghost: while conventional ghosts may have arbitrary large 
frequencies, and essentially, are introduced only to cancel unphysical polarisations of gauge 
fields for arbitrary large oj, typical frequencies of the Veneziano ghosts are order of the 
horizon scale, uj ~ H. 

The most immediate consequence of the interaction term (3.6) is that the magnetic field 
which will be generated this coupling will have a typical Fourier mode k EM ~ k of the same 
order of magnitude of the ghost mode. On the other hand, gravity can only lend energies 
of order u)k — k ~ H, higher frequencies being exponentially suppressed. Consequently, a 
typical EM mode will be around oj em ~ k E — H. Moreover, this interaction is active at 
every epoch, and we shall see in the forthcoming section how (3.7) naturally explains the 
/xG intensity apparently observed in galaxies and clusters, and maybe needed from the very 
beginning of structure formation. To conclude, let us notice once again that everything we 
have been discussing so far is part of the SM, and the coupling constants are all known. 
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However, all our conclusions will apply to most (pseudo) scalar models of dark energy, as 
long as they are augmented with a coupling with the same structure as (3.6); in this sense 
the results we will discuss in the upcoming section are general, and widely applicable to 
dark energy theories. 

4. Dark Energy and Large Scale Magnetic Field, or Two Sides of the Same 
Coin 

The Lagrangian density (3.7) explicitly contains a coupling between the ghost and the P- 
odd operator E ■ B made of electromagnetic fields. The structure of this term is identical 
to the textbook example describing the anomalous ir° — > 77 decay. Let us explore the 
consequences of this interaction in our specific circumstances. While the structure of the 
ghost-photon interaction is a photocopy of that of the pion-photon one, there is a funda- 
mental difference between the two: it is a massive physical particle which can decay to 
two photons, whereas the ghost is massless and unphysical, and can not decay into two 
photons. The ghost field in an expanding universe should be treated as the large correlated 
classical field which emerges from a non-zero expectation value ("H p hys|(02 — 0i)|%phys) 7^ 
as explained in [38]. The Fourier expansion of these classical fields fa and is saturated 
by very low frequencies uj^ ~ k ~ H , while higher frequency modes Uk S> H are strongly 
suppressed as a result of the relative suppression of the so-called Bogolubov coefficients [38] . 

4.1 Fast equilibration: first estimates of B 

For future convenience we introduce the dimensionless coupling constant which appears in 
our basic expression (3.7) 

(4.1) 

In nature /3 <C 1, but we still want to study the rate of energy transfer from dark 
energy (which is represented by the ghost field) to the electromagnetic energy as a function 
of j3. Hence, we treat (3 as a free parameter in this section. If the fields are interacting 
sufficiently strongly (/? > 1) and the potential minimum is reached sufficiently rapidly, 
one can estimate the expectation value (E ■ B) in terms of the external background field 
(4>2 — (/>i), which is treated as a source. By following this procedure we arrive to 

(E ■ B) = ^N f m q \(qq)\(^^) * ^ 3 QCB H , » 1 . (4.2) 

As expected, in Minkowski space there will be no generation of EM field as (02 — (f>i) = 
according to (3.4), (3.8). In an expanding universe {(f>2 — 0i) is proportional to the deviation 
from Minkowski space, and is expected to be around H. In this case the magnetic field is 
produced as an outcome of the energy flow from the ghost to the EM field. At the time 
at which each field mode is born, our equations are symmetric under the permutation of 
E and B. Therefore, one could estimate the absolute value of \B\ as 3 , (E ■ B) ~ (B 2 ). 

3 The exact statement is [{E 2 } + {B 2 )] > 2{E ■ B) but the electric field will be screened soon after it is 
generated, and it can then be neglected in the evolution. 
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Of course, the evolution of the electric field and magnetic fields are drastically different 
as electric charges do exist in nature while magnetic ones do not. However, we expect 
that (E ■ B) ~ {B 2 ) is a reasonable approximation for the absolute value of \B\ at its 
birth. Consequently, a simple estimate of the intensity of the magnetic field based on the 
assumption that the system (ghost + EM) can reach its minimum energy configuration 
sufficiently quickly (fast equilibration) can be presented as 

(B 2 )^(E.B)^^^^^, (4.3) 

where pde — Aqcd-^ in our framework. For large coupling constant /3 this estimate is 
totally justified. Indeed, eq. (4.2) corresponds to the minimum of the potential energy 
of the system 4 , and it becomes a precise statement for /3 > 1 when the equilibration is 
achieved faster than a Hubble time 1/H. 

However, in nature ^< 1. Let us see what is happening with our formula (4.3) when 
we start to decrease (3. This equation tells us that we generate magnetic energy which 
(parametrically!) starts to exceed the energy of the source when /3 ~ 1. Clearly this can 
not be true, and the loophole in the above line of arguments lies in the assumption of 
equilibrium. The minimisation procedure described above works only when the reactions 
involved in transferring energy from one source (the ghost) to a recipient (the electromag- 
netic field) are efficient enough to be in equilibrium. The lowest energy of the system 
simply can not be achieved when the coupling between the two components is weak. 

Thus, in a slowly expanding universe one needs to have a handle on the rate of energy 
transfer at each epoch, and compare it with the Hubble time, just as it is typically done 
in studying the thermal history of particle species in the early universe. The difference 
here is that this rate depends itself on the rate of expansion H due to the fact that the 
effective coupling constant is proportional to (fa — 4>i) — H as formula (4.2) states. As 
we will see shortly, in one Hubble time only a very small fraction of the ghost's potential 
energy is transferred to the magnetic field at f3 <C 1, thereby mining the foundations of the 
minimisation procedure. Yet, there is an important lesson to be learnt from this discussion: 
the dynamics described by the Lagrangian (3.7) does lead to energy transfer from the 
ghost field (which is the DE in our model) into the magnetic field, though numerically 
formula (4.3) can not be trusted for the physically relevant case /3 <C 1. 

4.2 The time scales 

As we mentioned above, in order to understand the dynamical effects of the coupling 
between the ghost and electromagnetism one should look for the rate at which energy is 

4 In particular, a similar procedure of minimisation of the effective potential allows to compute the exact 
vacuum expectation value for the gluon field (f^G^„G'" /a ) = -^m q {qq) which is analogous to (4.2) when 
the colour gluon fields G M " a replace the EM fields, and the so-called 9 parameter of QCD replaces the 
expectation value of the ghost field {(f>2 — 4>i). The equilibrium for (§jGji„G'"' a } is obviously achieved very 
quickly and the relation (fjG^G""" 1 ) = -^m q {qq) becomes exact at the minimum of the potential. In 
fact, one can differentiate this relation with respect to 6 one more time to arrive to well-known exact Ward 
Identity for the topological susceptibility, i J dx(^G%G^ a (x), G%& va {Q)) = ^f(qq), see, e.g., [60] 
and references therein. 
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transferred while the universe expands. In order to do so let us write down the Hamiltonian 
for the system (3.7) as 



H 



(4.4) 



+ 



a 
2^ 



X E- B - N f m q I (qq) | cos ' 



/n 



As we discussed earlier in the text, neither the Lagrangian (3.7), nor the Hamiltonian (4.4) 
contains any small coupling constant in their definitions as all the parameters describing the 
system are known SM parameters. All small effects are proportional to H/Aqqd ~ 10 -41 
and are brought about only at the level of the "renormalisation procedure", i.e., when 
the subtraction is explicitly performed. If we are computing the energy related to the 
magnetic field in Minkowski space, this identically vanishes (4.2) as a consequence of the 
conditions (3.4). 

The time derivative dH / dt governs the efficiency of the energy flow between ghost and 
electromagnetic fields. Explicitly, dH/dt = implies that 



_d 
dt 



1 



B 2 + E 2 ) 




E-B-/3 



fr> 



_d 
dt L 



EB 



dt 



-H 



,(4.5) 



where dH g host/dt essentially describes the dynamics of the dark energy component, and 
can be neglected in this simple evaluation. 

We want to gain some intuition about the rate of the energy transfer by considering the 
case of small (3 when we know that the minimisation procedure is not justified. Assuming 
for simplicity that the Hubble parameter is a constant and the rate of energy transfer is 
also a constant, one can estimate the typical time (relaxation period) to which is required 
for the system to reach its equilibrium value. Indeed, the left hand side of eq. (4.5) can be 
approximated as 



l.h.s. of eq. (4.5) 



d 
dt 



{B 2 + E 2 



PEN! 
TO 



At the same time, the right hand side of (4.5) is 



r.h.s. of eq. (4.5) 



h x ,92 



h > PEM 



fv 



(4.6) 



(4.7) 



where we have substituted (E ■ B) ~ pem and its time derivative with a 1/tq. The last 
step we need to do is simply to plug in the values (<f>2 — <t>i) / frf — {4>2 — 4>i) — H (non-zero 
solely due to the expanding gravitational background). A more precise numerical estimate 
for {<p2 — 4>i) will be given below. Comparing (4.6) with (4.7) one can immediately infer 
that 



to 



(3H ' 



(4.J 



which is the main result of this subsection. 
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4.3 Slow equilibration: more realistic values for B 

Now we want to see what happens with our estimate (4.3) when f3 decreases, as it corre- 
sponds to the physical value (4.8). From what we have said above, it is clear that the time 
scale which is required to attain equilibrium will be order of the Hubble time tq ~ ^-/H 
for f3 ~ 1 when the minimisation formulae can still be marginally trusted. In this case, 
a finite portion of DE of order one can be transferred into the magnetic energy during 
a Hubble time. When (3 decreases even further to reach the physically interesting region 
(3 <C 1 the formula (4.8) suggests that in order to reach equilibrium one needs a time which 
is /3" 1 ^> 1 times larger than the Hubble time 1/H, which really makes no sense. The 
appropriate interpretation in this case is that equilibrium will be never achieved with such 
small coupling constant during one Hubble time. Instead, a very small portion (3 of the 
available energy can be at most injected into the magnetic field within the same Hubble 
time. 

This is however not the end of the story for small /3 <C 1. The point is that, for 
(3 ~ 1, each event leading to equilibration transfers an amount of energy of order one. 
For small (3 <C 1 this argument does not hold anymore: not only the equilibration time 
increases as a result of fewer events per unit time, see eq. (4.7), also a smaller portion of 
the available energy will be transferred to the EM field per collision. It follows that a more 
realistic estimate shall account for an additional suppression factor f3 in comparison with 
case (3 ~ 1, when formula (4.3) can be still marginally trusted. 

With this interpretation in mind we arrive at our final estimate for the magnetic energy 
that has flowed from the ghost field 



where we inserted the extra small parameter (3 2 accounting for the two suppression factors 
mentioned above, accounting for the fact that only a small fraction of the total available 
energy will effectively exchanged during one Hubble time. 

One should remark here that the f3 2 suppression which enters in (4.9) can be intuitively 
understood by considering a system of particles with magnetic moment. In this case, as 
is known, the interaction between the magnetic moment and an external magnetic field 
is proportional to the coupling constant e. However, the contribution to the energy due 
to the induced magnetic moment is proportional to e 2 . In different words, the magnetic 
susceptibility of the system goes as e 2 rather than linearly with e, in spite of the fact that 
the strength of the interaction itself is oc e. Our analysis of the induced EM field follows the 
footprints of the analogous induced magnetic argument, where now the role of the coupling 
constant e is played by the parameter (3, and where the source of the energy is pde rather 
than the external magnetic field. Finally, in our expressions above we assumed that the 
Hubble constant does not depend on time, so that the transfer rate is time-independent 
itself. 

A few comments are in order before we proceed with the numerical evaluations. First 
of all, eq. (4.9) has the following parametrical dependence on fine structure constant a and 




(4.9) 
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Hubble constant H: 

Pem ot (B 2 ) ~ (^) 2 HA 3 QCD , where p DE ~ if • (4-10) 

In this form it is easy to interpret the appearance of each parameter. Dark energy is a 
QCD effect related to the mismatch between vacuum energies in infinite Minkowski and 
compact or curved spaces, and it is conceivably proportional to the rate of the expansion 
of the universe. As is known, the EM field does not interact directly with gravity, however, 
it does interact (through a standard triangle anomaly) with a (pseudo) scalar ghost field 
which gives rise to a correlated external source for the EM field in an expanding universe. 
This explains the extra parameter (a/2ir) 2 in (4.10). There are other numerical factors of 
order one which would also appear in (4.10), which for now are glossed over to simplify 
the interpretation. However, we believe that our formula (4.10) has a well understood and 
physically motivated behavior in terms of the relevant physical parameters a, if, Aqcd- 

We are now in the position to make a more quantitative estimate for the magnetic 
field as the combination which enters expression (4.9) and which describes the dynamics 
of the ghost field, that is, c ~ (fo — <j>i ) / (Hf v * ) , had been previously discussed in [38] 
and it essentially corresponds to the initial field velocity in the classical potential. An 
uncertainty in value of this constant, c, is not related to any physics beyond the SM, but 
is simply determined by the initial conditions when the dynamics of 02 and (f>i are treated 
classically 5 . For our order of magnitude estimate the specific value of this parameter is 
unimportant, and, following [38], we will fix it at the value of c ~ 10~ 3 . 

With these remarks in mind we can finally work out some numerology. The typical 
value for the energy density of the magnetic field which is generated by the Veneziano 
ghost, which in turn is responsible for the cosmological dark energy, is 

b 2 ^^N^Ql) 2 -c-pde, c-ht 3 . (4.ii) 

From this expression it straightforwardly follows that during the last Hubble time of life 
of the universe, an OiHQ 1 ) correlated magnetic field is born with intensity 

B ^ ^ v / c-(2.3-10" 3 eV) 2 ~nG, 1G = 1.95 • 10" 2 (eV) 2 . (4.12) 

As we mentioned previously, formula (4.12) should be treated as an order of magnitude 
estimate due to a number of numerical factors which have been consistently neglected while 
arriving at our final analytical result (4.9) and its numerical expression (4.12). We should 

5 The classical treatment of the system is by far the widest chosen approach in dealing with dark energy 
matters. However, while in other, perhaps more familiar, cosmological models such as inflation, the passage 
from quantum to classical is justified a posteriori (see for instance the discussion in [10], section 7.4.7, 
and reference to original works therein), in coping with our quantum fields we do not expect such a "little 
miracle" to happen. In other words, the quantum nature of our fields which appears in their non-trivial 
dependence on the gravitational background as well as on the global properties of the manifold, is brought in 
at the level of the renormalisation procedure, and is therefore not describable in a purely classical approach. 
Nevertheless, as long as we are interested in order of magnitude estimates, it will suffice to confine ourselves 
within the boundaries of such classical framework. 
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emphasise that the most important outcome here is not the precise numerical value found 
in (4.12), which is likely to change due to the very complicated evolution of the field in its 
most recent history (z ~ 1) well after its formation. Rather, the main result of this paper 
is the observation that due to the coupling with DE the EM field will be correlated on 
enormous scales of order 1/H. This distinguishes our mechanism from everything which 
has been suggested previously. Another important qualitative consequence is the prediction 
that parity will be locally violated on the same scales for which the field does not change, 
that is, 1/H as a result of the pseudoscalar nature of the DE field. This distinct and 
unambigious prediction can be tested in the CMB sky, where apparently parity violation 
on such scales indeed has been observed, see [61, 62] (we will comment further on this 
aspect below). Finally, our field is highly helical with typical wavelengths A& ~ 1/H; 
therefore, the induced B is expected to flip sign on the same scale. We notice in passing 
that a nG intensity would account for the observed galactic field if it were frozen in the 
pre-galactic plasma. In spite of its attractiveness however, this possibility should be taken 
with a grain of salt; furthermore, the analysis of the evolution of such fields is beyond the 
scope of the present work. 

This mechanism operates all the time, and there will be an uninterrupted flow of 
magnetic fields produced at different correlation lengths (proportional to the Hubble pa- 
rameter). Such fields however, being proportional to the vacuum energy component, will 
be of significantly relative weaker amplitude compared to the other components of the 
universe, following the twin evolution of the dark energy. Nevertheless, they could still 
behave as seeds for plasma mechanisms to process them, the outcome of which is beyond 
the scope of this paper and shall not be pursued any further. 

4.4 Some applications 

As we have already mentioned, interactions of the form (3.6) have been detailedly exploited 
and analysed in the literature, initially in the context of the anomalous axion-photon 
interaction, and then extended to general pseudoscalar interactions. What makes this work 
depart from all the paper referenced in the bibliography is that we are dealing entirely with 
SM physics, including the interaction (3.6). We should remark here that this interaction 
does not violate P and CP invariance on the fundamental level, similarly to the ir° — > 27 
decay. However, on small scales, one can interpret the interaction (3.6) as a P and CP 
violating coupling similar to the 6 term in QCD. Such a violation occurs only locally, on 
the correlation scales Xk ~ 1/H, while globally one should expect P and CP conservation 
according to the symmetry of the fundamental interactions (3.7). 

In this paper we have worked out the first application of such interaction in our model, 
that is, the generation of cosmological magnetic fields. It is easy to see, however, that 
there are several possible effects in addition to what has just been computed, especially 
in connection with possible observables which would be able to confirm or falsify our 
framework. We will limit ourselves to a mention for some of these possibilities, the specific 
details of which are beyond the scope of the present work and will be addressed in future 
publications. 
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• Early universe large scale magnetic fields have a sound impact on the mechanisms of 
structure formation (especially an active source of nanogauss intensity); there exists 
a vast literature on the subject, although this point is still somewhat overlooked: 
we refer to the comprehensive reviews [1] and [2] and to section 2.4 where some 
interesting consequences of primordial magnetism in the formation of early structures 
are investigated. 

• The presence of a pre-decoupling magnetic field is able to leave significant imprints 
in the cosmic microwave background radiation (CMB). The literature on this specific 
topic is particularly copious, and we shall mention only a few specific papers. The 
effects of magnetic helicity (particularly relevant for us since in our case the magnetic 
field is highly helical) have been discussed in [63, 64]. Primordial magnetism is also 
able to de-polarise the CMB (which acquired a small degree of linear polarisation 
due to the finite thickness of the last scattering surface), see [65, 66, 67, 68, 69]. Fi- 
nally, an impressive and complete analysis of the distortions caused by pre-decoupling 
magnetic fields (stochastically distributed) to the CMB acoustic peaks has been re- 
cently undertaken in the series of papers [70, 71, 72, 73, 74]. The novel aspect of our 
proposal compared to the models adopted in all the aforementioned efforts is that 
our source is active at all times, has a highly non-trivial redshift dependence, and is 
capable of producing otherwise forbidden by parity cross-correlation spectra such as 
those known as EB and TB on top of the distortions predicted for the usual TT and 
EE (and BB) autocorrelations and the parity even TE cross correlation. Anomalous 
parity violation in the low multipoles of the CMB has been discussed in [61, 62]. 

• In addition to the pervading magnetic fields generated through the Veneziano ghost, 
one could also look for other signatures of the parity- violating 6 coupling (3.6); in 
the early expansion history this interaction is able to trigger a sizeable amount of 
birefringence of the CMB [75, 76, 77] (the peculiar effects of the parity violation have 
been considered in [78, 79, 80, 81]); at late times such coupling can be sought after 
in its impact to the travelling light signals from the CMB, or any emitting source 
in the universe, to us, the observers, see again [49, 52, 54] and [87, 88, 89]. The 
latter is of prominent relevance due to the recent claims of a preferred orientation of 
the polarisation angles of quasars in optical wavelengths [82, 83, 84] which seems to 
disappear at radio frequencies [85] (see also [86] for an analysis of the possible biases 
affecting these results): a possible explanation is again to be found in the combined 
effect of a background magnetic field and an anomalous coupling of (3.6) fame which 
mixes light and the pseudoscalar field [87]. 

5. Results and conclusion 

In this work we have proposed that the origin of the observed cosmological magnetic fields, 
correlated over scales that stretch from galaxies to superclusters, is tied to the existence 
and evolution of the vacuum energy component of the universe. 

6 P and CP violation may occur locally, not globally, as explained above. 
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Here we have shown how the very same field which provides the dark energy in our 
framework [38], together with its Gupta-Bleuler partner, possesses an anomalous coupling 
(the ordinary triangle anomaly) to the magnetic field and can trigger the generation of 
magnetic fields of the observed intensity with correlation lengths of order of the Hubble 
parameter today. This interaction is present at all times starting from the QCD phase 
transition, and is continuously producing magnetic fields, whose strength is consequently 
linked to that of the vacuum energy, times a dimensionless coefficient of about (a/ Air) 2 
due to the weakness of the electromagnetic interactions. This combination, a 2 / 0D£/167r 2 , 
where in our set up pde — HAq CD , is precisely the magnetic field strength (squared) one 
needs to account for the observational evidence. 

The list of intricacies most large-scale magnetic field models have to face does not clash 
with the scheme we have outlined. Correlation lengths up to the size of the universe at 
each epoch are a direct consequence of the properties of the Veneziano ghost. Magnetic 
fields will be generated at all times, even at higher redshifts, though their magnitudes will 
experience a standard l/a 2 (t) suppression due to the expansion of the universe. There is 
no need for seeds to feed a dynamo with, as the largest scales are generated last, and with 
already the correct strength, without the risk of overwhelming BBN with gravitational 
waves. 

Distinctive predictions and unique features of our proposal are the generation of the 
magnetic field with magnitude in the nG range at all scales up to the Hubble radius today, 
and its grounds to be found entirely within familiar SM physics, without the need for any 
new fields, unconventional couplings to gravity or any modification of gravity itself. This 
mechanism is also unique in predicting parity violation on the largest scales 1/H, which 
apparently is already supported by CMB observations [61, 62]. If these predictions were 
to be confirmed by future PLANCK observations, it would be a strong (but still implicit) 
hint towards the QCD nature of dark energy, together with its fundamental kinship with 
large scale magnetic fields, in the way detailed in the present work. 

Our final comment concerns the scales involved in the model. Our proposal leads to 
the order of magnitude estimate B ~ ^ yj H Aq CD ~ nG, which is approximately the value 
that, by simple adiabatic compression, could explain the field observed at all scales, from 
galaxies to superclusters. It is already the second time that this "accidental coincidence" 
happens, the first one being the dark energy density itself, written as pde — H^qcd ~ 
(10~ 3 eV) 4 . We consider this "coincidence" as encouraging support for the entire framework. 
The SM coupling of the EM field with the DE field leads to the observed pG galactic 
magnetic field (assuming it were frozen in the pre-galactic plasma). It is difficult to invent 
another scheme when these relations hold, see, e.g., [90]. This "coincidental" relation 
should be added to our "Fine-tuning without fine-tuning" section from [38] where we clarify 
how the intimidating list of fine-tuning issues which always plagues dark energy models, 
possesses simple explanations without the need to introduce new fields, which come with 
new interactions, new coupling constants, and new symmetries. Instead, all our results 
are based on the paradigm according to which the "physical dark energy" is related to 
the deviation of the vacuum energy from infinite Minkowski space similarly to the Casimir 
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effect, whife afl SM fields couple to the dark energy field (the Veneziano ghost) in the well 
understood way dictated univocally by the standard model of particle physics. 
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